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ABSTRACT 



The mid-infrared (MIR) spectra of dense photodissociation regions (PDRs) 
are typically dominated by emission from polycyclic aromatic hydrocarbons 
(PAHs) and the lowest pure rotational states of molecular hydrogen (H 2 ); two 
species which are probes of the physical properties of gas and dust in intense UV 
radiation fields. We utilize the high angular resolution of the Infrared Spectro- 
graph on the Spitzer Space Telescope to construct spectral maps of the PAH and 
H 2 features for three of the best studied PDRs in the galaxy, NGC 7023, NGC 
2023 and IC 63. We present spatially resolved maps of the physical properties, 
including the H 2 ortho-to-para ratio, temperature, and G /nn- We also present 
evidence for PAH dehydrogenation, which may support theories of H 2 formation 
on PAH surfaces, and a detection of preferential self-shielding of ortho-H 2 . All 
PDRs studied exhibit average temperatures of ~ 500 - 800 K, warm H 2 column 
densities of ~ 10 20 cm" 2 , G Q /n H ~ 0.1 - 0.8, and ortho-to-para ratios of ~ 1.8. 
We find that while the average of each of these properties is consistent with 
previous single value measurements of these PDRs, when available, the addition 
of spatial resolution yields a diversity of values with gas temperatures as high 
as 1500 K, column densities spanning ~ 2 orders of magnitude, and extreme 
ortho-to-para ratios (Rop) of Rop < 1 and Rop > 3. 

Subject headings: ISM:clouds — ISM: lines and bands — ISM: photon-dominated 
region — ISM: individual (NGC 7023, NGC 2023, IC 63) — ISM: molecules — infrared: 
ISM 
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Introduction 



Photodissociation regions (PDRs) around hot stars are the boundaries between the 
ionized HII region and the surrounding molecular clouds. PDRs consist of layers of atomic 
and molecular emission with a chemical profile dictated by the ratio G /ns, where G Q 
is the strength of the local UV radiation field (912 - 2000 A) in units of the Habing 



flux (1.6 x 10 3 erg cm 2 s 1 



( iHollenbach fc Tielens 



Habingi (119681 )) and Uh is the local hydrogen density 



19971 ). The typical mid-infrared (MIR; 5-35 /xm) spectrum of a PDR 
is dominated by the stretching and bending modes of polycyclic aromatic hydrocarbons 
(PAHs) excited by the absorption of i ndividual UV photon s and by the lowest pure 



rotational lines of molecular hydrogen (ILeger fe Puget 



19841 ) . H 2 is formed on the surface 



of dust grains and is destroyed by photodisso ciation within PDR s. Its formation rate has 



been found to be related to PAH abundance (IHabart et al. 



20041 ). suggesting that there 



may be an observational relationship between these two species. Both PA Hs and H? have 



been shown to be valuable pro 



1994 



Draine fc Bertold 



1996 



jes of the physical conditions within PDRs fjBakes Tielens 



Hollenbach fc Tielens 



1997 



Berne et al. 



2009^ 



PAH emission in the MIR is dominated by the five broad emission features centered at 
6.2 /im, 7.7 (Urn, 8.6 /zm, 11.3 /mi and 12.7 /mi. These, along with features at 3.3 /im and 
17 /im, are also referred to as the Unidentified Infrared (UIR) features and the Aromatic 
Infrared Bands (AIBs). Studies have suggested that each of the major PAH bands is 
actually a blend of emission from PAHs of various sizes, surface structures, and charges, 
therefore the relative shapes, centroids and flux ratios with respe ct to the other bands 



vary depending on the physical properties of the observed PAHs (IBakes fc Tielens 



vary depending on t 


Peeters et al. 


2002 


• - 



1994 



van Diedenhoven et al. 



20041 ) . Laboratory measurements have shown 



that the features attributed to carbon-hydrogen (C-H) bond bending, namely the 3.3 /im, 
11.3 /tm and 12.7/im bands, are dominant for a population of neutral PAHs, however the 
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carbon-carbon ( C-C) stretching modes, from 6-9 nm, are the dominant e mission from PAH 



cations (PAH + ) ( lAllamandola et al 



1999 



Hudgins &: Allamandola 



19991 ). In recent years 



the ratio of C-C emission to C-H emission, most notably the 6.2 fim or 7.7 u m to 11.3 Aim 
bands, has become a meth od of estimating the charge state of a PAH cloud ( jJoblin et al. 



1996 



Galliano et al. 



20081 ) . A smaller number of studies have shown that the ratio of the 



12.7 (Ltm emission, which arises from duo and trio C-H out-of-plane bending modes, to the 
'. 1.3 um singlet C-H out - of-plane bending mode is a pro be of the PAH surface structure 



( IDulev fc Williams 



1981 



Hony et al. 



2001 



Tielens 



20081 ) . For smaller PAHs (with fewer 



than 50 carbon atoms), I12. 7 /hi, 3 m ay be an indicator of dehydrogenation, or the breakdown 
of duo and trio C-H bonds into solo bonds. 

In this paper, we present spectral line maps of the H 2 and PAH emission measured by 
the Spitzer Space Telescope from four of the best studied PDRs; NGC 7023, IC 63 and two 
separate regions in NGC 2023 (Figure s [D [21 [3]). All three PDRs are known to exhibit UV 



fluorescence of H2 ([France et al. 



(120091 ) and re ferences there i n). W e utilize the high angular 



resolution of the Infrared Spectrograph (IRS; iHouck et al.l ( 120041 )) onboard Spitzer to 



investigate the spatial variation of the temperature, column density, and ortho-to-para ratio 
of the lowest pure rotational states of H2, the ratio G /nn determined from the ionization 
state of the PAHs, and we address the possibility of dehydrogenation of PAHs across the 
observed regions. We present these maps and the methods used to generate them in §3. In 
§4 we present a summary of the new information revealed by the mapping technique for each 
PDR, and address the agreement between our measurements and previous observations. In 
§5 we address and similarities or differences between the PDRs and summarize our findings. 
A summary of the basic properties of the regions is presented in Table [TJ 
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2. Observations and Data Reduction 

2.1. Data Cube Assembly 

Spitzer data products were downloaded from the Spitzer Science Center (SSC) archive 
for NGC 7023, IC 63, and two locations in NGC 2023: NGC 2023 North (NGC 2023N), 
located ~ 160" north of the central star HD 37903, and NGC 2023 South (NGC 2023S), 
located ~ 78" south of HD 37903. IRS short wavelength low resolution (R ~ 60 - 127, SL) 5 
- 14 /im Basic Calibrated Data (BCD) spectra were obtained for each nebula. For a portion 
of NGC 7023 and for all of NGC 2023N we also obtained long wavelength low resolution (R 
~ 57 - 126, LL) 14 - 38 tan spectra. The se spectra were assembled into data cubes using 



the CUBISM software flSmith et al. 



20071) 



Spitzer BCD products require background and bad pixel subtraction. In IC 63, the 
outrigger exposures (the exposures of channel not centered on the target) were clear of 
nebulosity and therefore supplied a suitable background exposure. This was not the case for 
the extended PDRs NGC 7023 and NGC 2023. For NGC 2023S a dedicated sky exposure 
was taken and included in the BCD data package. To account for background emission in 
NGC 7023 and NGC 2023N, we searched the SSC archive for point source spectra of objects 
within 10° of declination of each of our targets taken within ± 3 days of our data. If the 
outrigger orders of these exposures were free of stars or nebulosity, then it was considered 
suitable for background subtraction. 

Bad pixels were removed using the CUBISM automated bad pixel routine as well as 
by interactive selection. Mapping observations with overlapping slits along the slit axis 
required editing of the default 'WAVSAMP' to eliminate contamination caused by the 
higher noise at the slit endpoints. The resulting 3-D data cubes (two spatial dimensions 
and one spectral, N x M x A) from each of the two IRS channels were merged into one by 
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averaging the regions of overlap between each of the spectral orders (the bonus orders were 
ignored). The spatial dimensions of the data cubes are ~ 1.8 x 1.8 "/pixel. A summary of 
the Spitzer AORs, program ID's and background program IDs for each of our datasets is 
listed in Table El 



2.2. Spectral Feature Fitting 



Each 1-D spectrum in the data cube was fit using the PAHFIT IDL routine (ISmith et al. 



20071 ). PAHFIT uses Drude profiles to fit the broad PAH features and Gaussians to fit the 
atomic and molecular lines, including the H 2 0-0 S(0) - S(7) rotational lines. The use of 
Drude profiles for the aromatic bands has been shown to recover ~2-6 times more flux in 



the broad wings o: 



( ISmith et al 



the 7.7 /zm and 8.6 /zm features than spline fitting or Gaussian profiles 



2006). The resulting integrated fluxes for each feature were then reassembled 



into maps. 

PAHFIT has a built-in routine to correct for extinction, including the silicate absorption 
features at 10 and 18 /im. To improve the spatial correlation of the extinction fitting, we 
re-binned our spectral maps into larger, 18 x 18" pixels (~ 5 x 5 IRS pixels). The resulting 
extinction values were extrapolated, smoothed and then applied to the original maps using 
the PAHFIT extinction curve. For the majority of our sample the MIR extinction is nearly 
negligible (the optical depth at 9.7 /mi, r 9 7 < 0.2), which is consisten t with the ISO deri ved 



global value inferred from observation of HD 200775 for NGC 7023 (IFuente et al. 



20001 ) 



and also with extrapolat ions of optical and near-IR measured extinctions into the MIR 



( iDraine fc Bertoldi 



20001) 



The data for NGC 2023N has an anomalous absorption feature between 9.3 and 10.7 
(Ltm with a position dependent centroid. The average full width at half maximum of the 
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feature is ~ 0.32 /zm and the average integrated flux is ~ 27.25 MJy /im sr _1 . The centroid 
varies between ~ 9.25 /im and 10.65 /im. We attempt to correct for by first determining 
the centroid and then scaling the original spectra over that region using a multi-gaussian 
fit and the average FWHM and integrated flux listed above. This correction introduces 
a systematic uncertainty to the H2 S(3) rotational line. Due to the importance of S(3) 
as the primary SL band ortho-rotational line detected in NGC 2023N, we utilize S(3) in 
our analysis despite the uncertainty. We estimate the error in the S(3) flux to vary from 
~ 0% near the edges of the observed region to up to 40% in the middle, however we see 
no position dependant artifacts in our temperature, column density, or H 2 ortho-to-para 
ratio maps (§3). Nevertheless, this ad hoc correction should be taken into account when 
considering the validity of our results for that PDR. 



3. Spectral Maps 

Spectral maps of the total H2 and PAH emission and the resulting line ratios and 
diagnostics presented in this section are included in Figures 0] - [7J Figures 0]- El (a) and (b) 
represent the total detected SL band PAH and H 2 emission respectively. Figures IH - [7] (c) 
and (d) represent the ratio of the 7.7 /im to 11.3 /im and 12.7 /im to 11.3 /im PAH bands 
respectively. Figures H] - [7J (e), (f), and (g) represent the calculated H 2 column densities, 
temperatures, and ortho-to-para ratios. Figures IH - [7J (h) show the calculated values of 
G /uh- In the following section we outline how each of these maps is constructed and 
discuss their physical interpretation^ 

For consistency, spectral analysis is only performed on data from the IRS SL channel, 



1 Maps of each individual H 2 rotation line, the major PAH features, and properties calcu- 
lated in this paper are available online at www.pha.jhu.edu/~flembri/PDRs. 
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except where noted. Therefore when we refer to PAH emission, we are neglecting the PAH 
features around 17 /mi, and when we refer to H 2 emission, we are neglecting S(0) and S(l). 
All maps have been oriented so that north is up. A summary of the average values for each 
of the properties calculated is presented in Table [3j These average values are not the same 
as the values obtained if the entire observed region was rebinned into a single pixel, however, 
as equal weight is given to each pixel regardless of the flux of the lines used to calculate 
that property. Therefore, the temperature of a pixel with an order of magnitude lower H 2 
flux than another is counted the same when computing the average, even though the pixel 
with the higher flux would dominate in a larger aperture. Single value measurements are 
already available for these PDRs and listed in Table [TJ 



3.1. PAH and H 2 Emission 

Individual spectral maps were constructed for the major Spitzer SL band PAH features 
(6.2, 7.7, 8.6, 11.3, and 12.7 /im) and H 2 0-0 S(J = 2-7) rotational emission lines using 
the extinction corrected line fluxes output by the PAHFIT routine. The signal-to-noise of 
the integrated line flux extracted from each 1-D spectra was required to be > 3 for the 
line to be included in its spectral map. Figures H] - [7] (a) and (b) represent the sum of the 
individual PAH and H 2 maps respectively and therefore are maps of the total detected PAH 
and H 2 flux in the 5-14 /im bandpass. 

The individual PAH features are strong, especially in NGC 7023 and NGC 2023, and 
Figures HI - [7] (a) represent the sum of the individual spectral maps with nearly complete 
coverage of the observed PDR in all bands. Figures H] - [7] (b), however, only represent a 
lower limit on the total SL band H 2 emission as the signal-to-noise of the higher J lines is 
often below our S/N > 3 threshold. Throughout this paper we neglect the H 2 S(6) (6.11 
/xm) line as it is heavily blended with the stronger 6.2 /im PAH feature and therefore cannot 
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be reliably extracted. 



3.2. PAH Band Ratios 



In recent years measur ements of th e PAH ionization state have become important 



diagnostic tools for PDRs. 



Berne et al 



(120071 ) successfully isolated PAH and PAH+ 
emission using blind signal separation techniques, showing that PAH + emission is most 
prevalent closer to the HII region where the radiation field is stronger. We investigate 
the degree of PAH ionization using the ratio of the 7.7 /im C-C stretching modes to the 
11.3 yum C-H out-of-plane bending mode. Laboratory studies have shown that ionized 
PAHs emit more in the C-C modes relative to the C-H bending modes than neutral PAHs, 
therefore Figures H] - [7] (c) are r e presentations of the ioniz ation state of the PAH molecules 



in the PDRs (jJoblin et al. 



1996 



Allamandola et al 



19991 ). The ratio I7.7//11.3 is expected 



to be ~ 1.3 for a 



( )Li fc Draine 



200l| ) . 



ully neutral column of PAHs and ~ 12 for a completely ionized column 



The maps presented in Figures H] - [7] (d) represent the surface structure of the PAHs 
given by the ratio /12. 7/^11. 3- PAHs can be stripped down to their carbon cores, or even 
dissoc iated via carbon fra gment loss, by absorbing UV photons with energies as low as 5 - 



8 eV ( ILe Page et al.ll2003l ). The dehydrogenation of a PAH decreases emission in the C-H 
singlet bands (centered at 11.3 fim) relative to the duo and triplet bands (centered at 12.7 
/im), however this relation is not fully understood, therefore we will focus only on what the 
trends in /12.7//11.3 reveal rather than the absolute values of the ratios. 



3.3. H 2 Temperature, Column Density, and Ortho-to-Para Ratio 



We can determine a gas temperature, T gas , using the pure rotational states of H 2 . 
These transitions have small critical densities and therefore can be assumed to be dominated 
by collisional de-excitation (see chap. 9.7.3 of Tielens (2005) for a review). Assuming 
thermodynamic equilibrium and a Maxwell-Boltzmann distribution, 



V gj+2 ' kT gas 
where Ej +2 is the energy and gj +2 is the statistical weight of the J+2 state. H 2 is a 
homonuclear molecule and hence has no dipole allowed pure rotational transitions, therefore 
we only detect quadrupole transitions with A J = 2. The statistical weight of the J state is, 



(2) 



<?J=0,2,4,6... — 2J + 1 
9J=1,3, 5,7... = (2 J + l)RoP 

where Rqp is the ortho-to-para ratio of the H 2 , w hich is 3 in a gas at temper ature greater 



than ~ 250 K in local thermal equilibrium (LTE) ([Sternberg fc Neufeldlll999l ). For H 2 out 
of LTE, the ortho-rotational lines (odd J states corresponding to parallel nuclear spins) will 
be offset from the para-rotational lines (even J states with anti-parallel spins), resulting in 
a "zig-zag" distribution. The column density, Nj +2 , of the J+2 state is determined by the 



Einstein A$(j) coefficient, 



47r/ s(J) 

7V J+2 - ~ A 7 77 , 1<3J 

where S(J) denotes the H 2 0-0 ( J + 2) — > J rotational transition and fsm is the measured 
flux density of the rotational emission line. This equation is valid only because the H 2 lines 
are dipole-forbidden and therefore optically thin. 
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In a uniform, single temperature gas with Rop = 3, the curve established by Equation 
[1] will be a straight line on a log-linear plot. The slope of this line is defined by the H 2 
rotational temperature. In a single temperature gas, the temperature, column density, and 
Rop can be determined by iterating Rop until the data can be fit by a single Nh 2 and T gas . 

PDRs consist of a multiphase gas that can be approximated by a low density, 



diffuse medium inte rspersed with high density clumps or filaments (IChokshi et al. 



Sellgren et al 



1988 



19921 ). The majority of the H2 is in the cool, low density phase, however the 
emission spectrum in the SL band is dominated by the higher density warm component, 
often making up < 1% of the total H 2 column (Figure [S]). For the purposes of this study, 
we assume a single temperature dominates the H 2 S(2) - S(5) rotational emission lines 
and fit a temperature, N H2 and Rop accordingly. To construct maps of these properties 
(Figures H] - [7] (e), (f) and (g) respectively), we require that at least three of the four H 2 
S(2) - S(5) rotational emission lines have a signal-to-noise ratio greater than 3 and that the 
line established by the fit falls within the la error bars of all three of the measurements. 
Due to the lack of LL coverage in our sample and the low statistical weight of the S(0) line, 
establishing the temperature, N H2 , and Rop of the cool H 2 is not possible except in NGC 
2023N and a limited portion of NGC 7023. 



3.3.1. H 2 S(4) Overestimation 

In NGC 7023 and NGC 2023, H 2 S(4) consistently falls above the fit defined by S(2), 
S(3) and S(5). We determined that this is due to the S(4) line, at 8.06 /zm, overlapping the 
broad 7.7 /im PAH feature. The H 2 features in both of these PDRs are weak compared to 
the PAH emission (see Figure H] - (a) and (b)) and the PAHFIT routine overestimates 
the S(4) flux in these two PDRs by ~ 10-50%. This error causes an overestimation of the 
H 2 temperature of ~ 30 K and N(H 2 ) of ~ 10%. For each PDR, we used our temperature 
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fit to the S(2), S(3) and S(5) lines to estimate the actual S(4) flux and then applied the 
average overestimation, relative to the PAH 7.7 fim flux, as a correction to the S(4) flux for 
the entire PDR. This correction was necessary as the S(5) line is weaker than the S(4) line 
in these PDRs and not always detected. 

The ortho-to-para ratio is much more sensitive to the S(4) overestimation than the 
temperature or column density. While our correction does improve the quality of the Rop 
fit, it also introduces an additional degree of uncertainty, therefore we do not use H2 S(4) to 
estimate ortho-to-para ratios. We utilize all S(2) - S(5) rotational lines to fit a temperature 
(which simultaneously fits an ortho-to-para ratio) and then we remove the S(4) line, set 
the temperature and column density, and determine the best fit ortho-to-para ratio to the 
remaining lines. This significantly reduced the scatter of the ortho-to-para maps (Figures H] 
~E] (g)). IC 63, which has lower total PAH flux relative to H2 emission (Figures [7] (a) and 
(b)), does not show evidence of H 2 S(4) overestimation. 



3.4. G /n H 

The gas temperature and the ionization state of PAHs can be used as a probe of the 
local density and UV field strength, G /uh- Using the parameter fit for t he 11.3 tan flu x 



ratio of ionized to neutral PAHs (^113/^113 ~ 0.6) and Equation [2] from 



Joblin et al. 



f ll996[ ). 



T ,\ _L K Z.A PAH+ 

J-7.7fim _ ^ 3 ^ PAH 

hl.^m \ 1 + 0.6^§^ 

where PAH + and PAH represent the total number of ionized and neutral PAHs respectively. 
This equation can be related to \jT gas G Q /n e by solving for the neutral fraction, f n = 
PAH / (PAH + +PAH° ), and setting it equal to the ionization parameter relation, f n = 



13 



(l+lo)' 1 ( Bakes & Tielenslll994l : lTielensll2005l ) 



7o 



3.5 x 10 _D A^- 



.lG n T 1 ' 2 



n, 



(5) 



7.2-0.6(/ 7 . 7 //n. 3 )/ c 



n, 



(6) 



We adopt an average number of carbons atoms per PAH of Nc = 50. This approximation 
breaks down near the total ionization limit of /7.7//11.3 — >■ 12, however we do not encounter 
J7.7//11.3 > 11 except in NGC 7023 in the inner nebular regions near HD 200775. We do not 
include data from that region in the following analysis. We estimate the electron density, 
n e by assuming that all free electrons come from singly ionized gaseous carbon, and that 



Sofia 



all ato mic carbon in the PAH emission region is ionized. Using the C/H ratio from 
(120041 ). the electron density is estimated to be n e ~ 1.6 x 10~ 4 ia H . Maps of G /nn are 
presented in Figures H] (h) - [7] (h) . 



Results 



4.1. NGC 7023 



In Figures H] (a) and (b) it is clear that the PAHs and H 2 do not share the same 
morphology. The PAHs extend ~ 20" closer to HD 200775 than the H 2 , which rapidly 
falls from its peak at the starwards edge of the ribbon feature. Figure [9] shows that 
the PAH emission peaks on average ~ 1.5" starwards of the H 2 . The separation of H 2 
rotational emission and PAHs has bee n shown before, most notably across the Orion Bar 



( jTielens et al. 



1993 



not been observed 



Kassis et al. 



Compiegne et al. 



2006), however in PDRs with weaker UV fields this has 



20071 ). 
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4.1.1. PAH Band Ratios 

Figure H] (c) shows that the PAH band ratio I7.7//11.3 is at a minimum near the bright 
ribbon feature and increases with decreasing project ed distance to HD 200775. This is 



consistent with the findings of 



Rapacioli et al 



( 120051 ) . who used ISO spectra to show that 



the relative contribution to the PAH spectrum of ionized PAHs increases closer to the 
central star. The geometry of the band ratio map does not match any of the other MIR 
observables and we note no apparent correlation between total PAH emission and ionization 
until I7.7/Z11.3 ~ 8. At this ratio, which from Equation [4] indicates an ionization fraction of 
/+ ~ 0.74, the total PAH flux begins to rapidly decline, potentially indicating destruction, 
falling by more than half over ~ 5" (Figure [H]). 

The dehydrogenation ratio /12.7//11.3 is highest both in the bright ribbon feature and 
close to HD 200775 (Figure |4] (d)). Beginning at the starwards edge of the dissociation 
front and extending ~ 10" towards the central star there is a trough where the ratio reaches 
a minimum. As with the ionization ratio, the geometry does not follow that of any other 
observable. 



4-1-2. H 2 Features 



The increase in H2 column density (Figure H] (e)) and flux at the bright ribbon feature 
is likely due to the onset of H 2 self-shielding, which is a rapid process. In the ribbon feature 
the Rqp of the H? is a nonequilibriu m value of ~ 1.83 ± 0.51, which is consistent with the 



2 reported in 



Fuente et al. 



(120001 ) . We note a spike in the ratio to values > 3 near the 



starward edge of the diss ociation front which is p ossibly the result of the greater ability 



of ortho-H2 to self-shield ( iDraine &: Bertoldi 



19961 ) . This measurement is not statistically 



significant, however, and is consistent with an Rqp = 3. 
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We observe a small but still significant separation of ortho- and para-H2 ~ 2", which 
is shown in Figure ITU1 Between 30-40" from HD 200775 there is an increase in Rop which 
occurs at the HI/H 2 transition. While the detection of Rop > 3 in NGC 7023 is not 
statistically significant, these profiles and the Rop spike suggest that we are observing the 
effects of preferential self-shielding of ortho-H2. If this interpretation is correct, it would be 
the first detection of the effects of preferential self-shielding of ortho-H 2 in a PDR. 



The bright ribbon of NGC 7 023 has been a target of numerous H 2 studies from t 



far-UV (France et al. 



2006 



and MIR f lFuente et al. 



lie 



20091 ) to the NIR flLemaire et al 



1996 



Takami et al. 



20001) 



20001 ) . however to our knowledge this work is the first to derive 



properties of H2 emission between the dissociation front and HD 200775, where the gas 



is subject to a harsher UV radiation field and higher G q 



1-0 S(l) and other narrow band imaging ( 
to detect any H 2 emission in this region 



Lemaire et al 



Fuente et al. 



Figure H] (h) ) . Prev ious H 2 



1996 



Takami et al. 



2000|) failed 



O2000|) reported a weak 0-0 S(l) 



detection with ISO, however they did not detect any other rotational lines at that location. 
We find that the IRS SL band H2 in this region is generally hotter (~ 1000K), at nearly an 
order of magnitude lower N H2 (N H2 ~ 1 x 10 19 cm -2 ), and close to an Rqp of 3. 



Previous MIR studies of NGC 7023 with ISO flFuente et al. 



2005|) and Spitzer flBerne et al 



2000 



Rapacioli et al. 



20071 ) have indicated a warm-H2 at T gas ~ 400 - 700K and 



Nh 2 ~ 10 20 cm -2 . These values are not matched by those listed in Tabled however the 
data is the table is skewed by giving equal weight to T gas and Njj 2 calculated for the H 2 
starwards of the bright emission ribbon. When the flux of each rotational line over the 
entire observing area is rebinned into a single value, which is dominated by flux from the 
ribbon, we calculate a T gas ~ 650K, Njj 2 ~ 6 x 10 19 cm -2 , and Rop ~ 1-83. These values 
are in line with previous work. 
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4.2. NGC 2023 North 

The PAH emission in NGC 2023N peaks ~ 10 " away from the peak of the total 
detected SL band H 2 emission. The separation is more pronounced between the PAHs and 
emission from the para-H 2 states (including S(0)) than the ortho-H 2 states (including S(l)). 
This can be seen by noting that the Rop peak in Figure [5] (g) is closer to the PAH emission 
peak in Figure [5] (a) than the H 2 peak in Figure [5] (b). This indicates that the direction of 
the advancing photodissociation front is moving from the lower right of the maps towards 
the upper left, even though HD 37903 is located closer to the bottom left. 



4.2.1. PAH Features 

The PAH ionization and dehydrogenation relations in Figures [5] (c) and (d) both reach 
a minimum near the PAH emission peak and increase with decreasing PAH emission. The 
low /7.7//11.3 values indicate that the PAHs are primarily neutral throughout the observed 
area. It is likely that the slight increase in PAH ionization to the north of the dissociation 
front in NGC 2023N is due to a decreasing n# and relatively flat G Q , therefore increasing 
the ratio G /nn as seen in Figure [5] (h). 



4.2.2. H 2 Features 



I 

NGC 2023N is also a well studie d PDR, with observations from the FUV (IBurgh et al. 



2002 



France et al. 



20091 ) to the NIR flHasegawa et al. 



1987 



Gatley et al. 



1987 



Takami et al. 



20001 ). Unlike NGC 7023, we calculate a nearly flat and featureless gas temperature and 
only a factor of ~ 2 drop in iV# 2 outside of the H 2 bright ribbon (Figures [5] (f) and (e)). 
When all of the SL band H 2 is summed, as in Figure [5] (b), the para-H 2 dominates the flux 



and N H2 . 
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The separation between ortho- and para-H2 is larger in NGC 2023N than in NGC 
7023. Figure [TT] shows a ~ 10" separation between the ortho- and para-H 2 peaks which 
exists not only in the SL band rotational lines (S(2) - S(5)), but also in the LL rotational 
lines S(0) and S(l). NGC 2023N is not oriented so that the projected vector towards the 
central star is perpendicular to the dissociation front, nevertheless Figures |5] (f), (g), and 
Figure HU show that Rop is nearly identical on either side of the Rop peak and T gas is 
flat throughout. This implies that collisional conversion has transitioned the Rop from a 
formation value to ~ 1.25 (Rop calculated using only the S(0) and S(l) lines for T gas = 
200 K is 0.67), and that the passing of the dissociation front is causing a temporary peak 
in Rop by preferential self-shielding of the ortho-H2 in both the cool, diffuse gas and the 
warm, dense clumps simultaneously. 



4.3. NGC 2023 South 



NGC 2023S contains a bright e mission ridge ~ 78 " south of H P 37903 which has been 



thoro ughly studied from the far-UV (IBurgh et al. 



20021 ) to the radio (ISteiman-Cameron et al. 



1997) and has been used as a template for modeling of dissociation fronts ( iDraine fc Bertoldi 



19961 ). The ridge area is considered to be denser (~ 10 5 cm~ 3 ) and have a slightly more 
thermal H 2 vibrational spectrum than its northern counterpart. The NGC 2023S IRS 
coverage extends over a significantly broader area than the bright emission ridge, however, 
extending up to ~ 15" from HD 37903 to as far as 100" away. 

The H2 and PAH emission in Figures [6] (a) and (b) both are at a maximum in the 
bright emission ridge centered at approximately [0,-20]. As with NGC 7023, the PAH 
emission remains nearer to its peak flux level at smaller projected distances to the central 
star than the H 2 . Figure [12] shows that the PAH emission peak precedes the H 2 emission 
peak in the bright ridge, as it did in NGC 2023N and NGC 7023, but by only ~ 0.6", which 
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is not statistically significant. Higher resolution measurements would be needed to confirm 
this separation. 

4.3.1. PAH Features 

The PAH ionization band ratio I7.7//11.3 in the bright emission ridge is at a similar 
average value (/7.7//11.3 ~ 3, /+ ~ 0.25) to the ratio in the ribbon of NGC 7023. To the 
north of the emission ridge the ratio increases to a maximum value of just less than 7, 
indicating an ionization fraction of / + ~ 0.65. In the mapped region closest to HD 37903, 
however, the PAH ionization level drops, possibly due to the complex geometry of NGC 
2023S resulting in foreground or background contributions to the PAH emission. 

The PAH dehydrogenation band ratio map (Figure [6] (d)) may be contaminated by the 
presence of a bright star at the bottom edge of the mapped area. The star and another 
fainter one to the north were removed to avoid contamination, however at 12.7 /im some 
contamination still exists. The concentric rings around the blacked out area in Figure [6] (d) 
suggest that this may be affecting the measurement. In the northern part of the map we 
observe values similar to NGC 7023 and NGC 2023N with slightly higher ratios nearer to 
the central star. 

4.3.2. H 2 Features 

The bright emission ridge coincides with the highest Njj 2 in NGC 2023S (Nh 2 ~ 1-87 ± 
0.30 x 10 20 cm -2 , Figure |6] (e)). Within this region the temperature averages ~ 625 ± 13K 
with an Rop ~ 1.87 ± 0.1. The spread on the values cited represent the standard deviation 
within the region, not the error on the measurements. We cite standard deviations here to 
show the relative homogeneity of of the observations within the regions discussed in this 
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section. Farther to the south of the ridge, Nh 2 drops off slightly to ~ 8.9 ± 3.5 x 10 19 
cm -2 , while T gas (592 ± 32K) and Rop (1.7 ± 0.2) remain at levels similar to those within 
the emission ridge. In the top left portion of the maps, nearest to HD 37903 and coincident 
with the area where we noted a drop in /7.7//11.3 in §4.3.1, Nh 2 , T gas and Rop are all 
consistent with the values in the southern portion of the map ( T gas ~ 575 ± 42K, Nh 2 ~ 
6.9 ± 1.6 x 10 19 cm -2 , Rqp ~ 1.63 ± 0.41). The similarity in these properties suggests 
that the gas dominating the emission in the upper left may be a foreground or background 
portion of NGC 2023S enveloping the central cavity. 

North-west of the dissociation front, Nh 2 drops by over an order of magnitude relative 
to the ridge feature {Nh 2 ~ 1.9 ± 1.4 x 10 19 cm -2 ). The corresponding increase in T gas 
to ~ 847 ± 325K suggests that the H 2 in this region of NGC 2023S is similar to the H 2 
starwards of the dissociation front in NGC 7023. The increased PAH ionization fraction and 
higher G /nn in Figure [6] (h) also follow the PAH behavior starwards of the dissociation 
front in NGC 7023. Rop, however, rather than increasing to the LTE Rop ~ 3, falls to 
extremely low values ~ 0.6 ± 0.22. This unusual observation will be discussed further in § 
5.3.1. 

4.4. IC 63 

IC 63 features the weakest PAH and H 2 flux in our sample and the lowest ratio of 
PAH flux to H 2 flux (Figure [7J (a) and (b)). We detect no statistically significant separation 
of the PAH and H 2 emission in IC 63 (Figure IT3|) . We also calculate the lowest average 
/7.7//11.3 ratio in our sample with an estimate ionization fraction of /+ ~ 0.082. This ratio 
does not reach the average values observed in the other PDRs anywhere in the observed 
region (Figure [7] (c)). It is worth noting that IC 63 is at a larger distance from the exciting 
star 7-Cas (~ 1.3 pc) than NGC 2023 (~ 0.33 pc to NGC 2023N) and NGC 7023 (~ 0.1 
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pc), therefore G D is lower due to geometric effects despite the harder radiation spectrum 
from the BO star. Figure [7] (d) shows that the dehydrogenation relation /12.7/A1.3 is mostly 
constant, however there is a decrease in the ratio near the starward edge of the PDR. We 
will investigate this further in §5.1.2. 

The temperature, column density, and G Q /n H all fluctuate only moderately (Figures [7] 
(e), (f) and (h)). The avera ge m easured TL S of ~ 625K is consistent with the temperature 



cited in 



Habaxtetal 



(120041 ) and 



Karr et al 



an extensive multiwavelength study in 



f|2005 Hr, 



Karr et al. 



- gas 



630K). IC 63 was the subject of 



(120051 ) in which ISO MIR data was use 



to compute temperatures and column densities, as well as to compare PAH band ratios. 
They reported a nearly flat PAH In %/I%.i band ratio, a probe of PAH ionization, which is 
in good agreement with the low variation of /7.7//11.3 shown in Figure [7J (c). They calculate 
an average In, 3/I&.2 ~ 2, however, which is twice the average that we calculate for this 
PDR (/11.3//6.2 ~ 0.93), and would indicate even less PAH ionization if correct. The Rop in 
Figure [7J (g) peaks at a nearly LTE value at the starwards edge of IC 63 before moderating 
down to ~ 1.5 at the back edge of our observing region. This implied trend in Rop can 
only be confirmed with a more extended coverage area, however it is consistent with an 
advancing photodissociation front. 



5. Discussion 

In this section, we address connections that can be made between the properties 
that we have observed in these PDRs. We address the limited understanding of the 
dehydrogenation relation and reference a proposed H2 formation scenario which may explain 
the observed band variations. We also consider whether our calculation of G /rin is a valid 
one given that the dominant H 2 emission in the IRS SL bandpass is from warm, dense 
gas clumps whereas the dominant PAH emission may not be. We conclude by analyzing 
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our ortho-to-para ratio maps and how observed trends fit with current understanding of 
advancing dissociation fronts. 



5.1. PAH Dehydrogenation 



5.1.1. NGC 7023 and NGC 2023 



In NGC 7023, NGC 2023N, and NGC 2023S, we observe an increase in the PAH 
dehydrogenation relation, 112.7/^11.3, beginning within the bright H2 ribbons and increasing 
slightly with greater distance from the central stars. This increase indicates that the PAHs 
are being stripped of hydrogen atoms despite corresponding drops in the ionization relation, 
/7.7//11.3, and G /uh (Figures HJ O [6] (c), (h)). In all three PDRs, /12.7//11.3 reaches a 
minimum, indicating the peak of PAH hydrogenation, just starwards of the H2 peak, where 
the H2 is dissociating. The /12. 7/^11. 3 relation then increases again nearest the central stars 
in NGC 7023 and NGC 2023S. This effect is illustrated in Figure H which shows the 
averaged profiles of /12.7//11.3 and N(H 2 ) with respect to HD 200775 in NGC 7023. The 
corresponding drop in total PAH emission near HD 200775 seen in Figure |9] may be evidence 
of the evaporating of PAH molecules due to a higher G /nu and ionization fraction. This 
would be consistent with previous results showing that PAH ionization precedes dissociation 



in the reflection/emission nebula IC 405 ([France et al. 



20071 ). 



The PAH dehydrogenation in the molecular layer may be due to a different process. 



H 2 formation on PA 



i surfaces has been proposed as a dominant formation process in high 



temperature PDRs ( iHabart et al. 



20041 ) . overtaking the traditional ISM formation process 



between physi s orbed H atoms on grain surfaces, which is effective only at low temperatures. 



Le Page et al 



( 120091 ) propose a formation mechanism where an ionized, fully hydrogenated 



PAH encounters an H atom and becomes overly hydrogenated (the gas phase hydrogen 
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atom becomes bound to the surface of the PAH) and then recombines with a free electron, 
which dislodges either the extra H atom or an H 2 molecule. In the H 2 ejection scenario, the 
neutral PAH product would be left with one fewer hydrogen atom than it had previously, 
and thus is dehydrogenated. The observed increase in /12.7//11.3 may indicate the presence 
of dehydrogenated PAHs due to an ongoing H 2 formation process. 

Neither PAH destruction nor H 2 formation fully explains the decrease in /12.7/J11.3 just 
starwards of the dissociation front, which is most easily seen in NGC 7023 (Figure H] (d)). 
This region has a higher PAH ionization fraction than in the molecul ar layer, which is an 



important catalyst for most H ? formation sc enarios on PAH surfaces (IBauschlicher 



Pauzat fc Ellinger 



2001 



Le Page et al. 



1998 : 



20031 ). therefore we would expect the formation 



rate to increase as more PAH cations become available. If the dehydrogenation ratio in 
the molecular layer is due to H 2 formation, then the lower I12.7/I11.3 starwards of the 
dissociation front would indicate a quenching of that process. 



5.1.2. PAH Destruction in IC 63 



We detect no significant variation of I12.7/I11.3 m the emission nebula IC 63 other than 
a minimum near the southern starwards edge (Figure [7] (d)). This region corresponds the 
highest Rop in IC 63, suggesting selective self-shielding of ortho-H 2 and therefore an HI/H 2 
transition layer. A dehydrogenation minimum near the leading edge of the HI/H 2 transition 

is consistent with observations in NGC 7023 and NGC 2023. IC 63 also exhibi ts FUV Y. 

1 



fluorescence indicative o f very little extinction between the PDR and 7-Cas (IWitt et al. 



1989 



France et al 



20051 ) and the PAH em ission region has 



an ionized atomic layer ~ 0.002 pc thick (IKarr et al. 



2005 



3een shown 


be ( 


France et al. 


2005) 



PAH emission does not continue into the HII region traced by Ha indicates that the PAH 
molecules are being dissociated near the starwards edge of the PDR. Figure [7] (d), however, 
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does not show any evidence of an increase in the dehydrogenation relation, suggesting 
that either that the dissociation is happening too quickly to be resolvable or that the 
dehydrogenation signature is too weak to detect. PAH molecules can be dissociated directly 
by the absorption of single high energy photons, however the required photon energy 
increases with the number of carbon atoms in the core, wit h a PAH of Nc = 50 requiring 



an EUV photon of ~ 24 eV (jSiebe nmorg en fc Krugel 



2010 



). Recent studies have shown 



20101 ). making single photon 



that 7-Cas is a soft x-ray emitter (ILopes de Oliveira et al. 
dissociation a possible scenario. 

The average value of the /12.7/A1.3 dehydrogenation ratio in IC 63 is similar to the 
average in the H 2 bright ribbons of NGC 2023 and NGC 7023 (~ 0.41), suggesting that 
the same processes that influence this ratio are occurring in all four PDRs. The ionization 
relation however (Figure [7] (c)), indicates that IC 63 is very close to the PAH neutral limit 
of /7.7//11.3 ~ 1.3, limiting the number of PAH cations available for H 2 formation. 



5.2. Reconciliation of G /nn 



To calculate G /nn we have used the PAH 7.7 jum and 11.3 /im flux in conjunction with 
the temperature of the warm H 2 component, wh ich represents ~ 1% of the total H 2 colu mn 



(Figure |8]). This method has been used before ( [Galliano et al 



2008 



Berne et al. 



20091 ). 



however it assumes that the PAH ionization relation is independent of the local hydrogen 
density. In our example, where we have used the warm H 2 temperature to calculate G /nn, 
we are assuming that we are calculating the ratios within those dense gas clumps. 

Radio observations of PDRs often use other molecules, such as CO, as a proxy for 



H 2 (IChokshi et al 



19881 ). These studies predict H 2 temperatures for the low density, cool 
H 2 component (n# ~ 10 3 cm' 3 , T gas ~ 100K), while MIR studies of the warm, dense H 2 
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clumps have indicated densities ~ 10 4 - 10 5 cm" 3 and temperatures ~ 600K. For these 
cases, the difference in temperature, \/T warm /T cool ~ 2.5, does not account for the order of 
magnitude or more difference in n# using the equations in §3.4. 

In order for the actual value calculated to be the correct G /nn within the warm, 
dense H 2 clumps that we are studying, the PAH flux must be dominated by emission from 
within those same high density clumps, thereby probing the same environment as the warm 
H2 temperature that we are using. If this is not the case, then either G D is larger within the 
high density clumps than in the low density, cool component, or we are improperly using 
the warm T gas and obtaining an incorrect measure of G /nn- We therefore caution that 
our maps of G /uh be used primarily as a measure of the variation of the ratio over the 
PDR, and not as a means of obtaining an absolute value for a specific portion of the PDR. 



5.3. H, Ortho-Para Ratio 



The formation ratio of ortho- to para-H 2 {Rop) i n a warm gas (T gas > 250 K) is 
expected to be 3. H 2 that forms at a cooler temperature will have a lower Rqp as more H 2 
is formed in the J = para-rotational state. Subsequent radiative heating of the gas allows 
J to change by only ± 2, thus the formation Rop is preserved in higher J rotational states. 
Radiative processes cannot convert ortho-H 2 to para-H 2 , however the total nuclear spin can 
be altered by collisions with H and H + atoms, whi ch moves Rqp towards the local thermal 



equilibrium value defined by the new temperature ( Sternberg Sz Neufel 



1999 ). 



The observed Rqp for a temperature, T gn , x , and density o: 



Neufeld et al 



' atom ic hydrogen, n HI) after 



(120061 ). Assuming a T, 



gas 



a time r can be predicted by Equation [2] from 
650K, a density of n# ~ 10 5 , uhi = 10 -3 n#, and an initial ortho-to-para ratio of Rqp 



0.75 (corresponding to a formation temperature ~ 50 - 75K), the observed Rqp after 5000 
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years at T gas should be ~ 2.15. This is within the la errors of the average Rop we measure 
for the ribb on of NGC 7023 (1 .86 ± 0.32), which has been exposed to HD 200775 for ~ 



5000 years ([Rogers et al 



1995|). It is also consistent with IC 63 (R OP = 1.99 ± 0.37) which 



has similar T gas and n#. 

A PDR with a moving dissociation front should exhibit a gradient in both the 
temperature and Rop as the front advances into the cool, low Rop gas, heats it, and begins 
the para — > ortho conversion. This process is most clearly observed in NGC 7023 in Figures 
H] (f) and (g), where the cool, outer regions of the PDR have a lower Rop than nearer to the 
central star, where Rop is closer to 3. It is also apparent in Figure [LTD where beyond the H 2 
emission peak R Q p falls to ~ 0.6, reflecting a low formation temperature, while T gas levels 
out, indicating recent heating. We also observe a spike in Rqp at the starwards edge of the 
dissociation front in NGC 2023N, however the measured values do not exceed 3. There is 
also the possibility that the Rop ~ 3 closest to HD 200775 in NGC 7023 is the result of 
rapid dissociation and formation processes and the relatively small amount of H 2 detected 
in that region is newly formed H2, not relic H2 which has transitioned to LTE. 

All four PDRs studied in this paper exhibit non-equilibrium ortho-to-para ratios which 
indicate that the H 2 formed at much lower temperatures prior to illumination by the central 
stars. The ortho-to-para ratios in IC 63 and NGC 7023 generally follow expectations for a 
recently heated PDR, however more coverage area in IC 63 is needed to confirm the implied 
trends in Figure [7] (g) . 



5.3.1. Extremely Low R OP m NGC 2023 South 



In NGC 2023S we calculate very low non-equilibrium ortho-to-para ratios within a 
pocket of warmer, low N(H 2 ) gas at a smaller projected distance to the central star than 
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the bright emission ridge (Figure |6] (g)). Ratios in this region, where the temperature 
occasionally exceeds 1000K, fall to < 1 with an average of 0.6 ± 0.22. At such high 
temperatures the para — > ortho conversion process should bring the gas to almost LTE 
within 1000 years an d to Rqp > 1 within a few hundred years for n# ~ 5 x 10 4 cm -3 



f lNeufeld et al. 



2006( 1 . We know of only two scenarios which could explain such extreme 
values. Either the gas has only been recently and rapidly heated, or the density of the 
gas is significantly lower than in the rest of the PDR, slowing the transition process. This 
region coincides with the highest G /nn in NGC 2023S (Figure [6] (h)), implying either a 



strong UV radiation field or a low density. S uch 



500K) have been reported in 



Neufeld et al. 



ow or tho-to-para ratios f or warm (T 



(J2006|) and 



Maret et al. 



)9) 



gas 



> 



(120091 ). however only 



for gas recently heated by a passing shock wave. This region merits further study to better 
understand and confirm these results. 



6. Conclusions 

In the photodissociation regions near hot stars, the MIR spectrum is dominated by 
emission from PAHs and the pure rotational lines of H2. We have shown that while these 
two species are independent of each other, they provide complementary diagnostic tools 
to probe the environment of the PDR. Utilizing the spatial and spectral capabilities of 
the Spitzer Space Telescope, we have created spectral maps of the ortho-to-para ratio, 
temperature, column density of H2, G /nn, PAH ionization, and a measure of the PAH 
surface structure. These diagnostics, combined with spectral mapping techniques, enable 
the study of PDR environments on angular scales of a few arcseconds. Such maps will prove 
valuable to future spectral mapping studies with Herschel, which does not have direct 
access to the bulk of the PDR gas mass (H 2 ) in its far-infrared and sub- mm bandpasses. 
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6.1. PAH Features 

Using the PAH emission in the MIR, we were able to establish an estimate of 
\/Tg as G /nii across the entire coverage region using the h.inm/h\.z^m C-C/C-H band ratio. 
There is evidence for PAH dehydrogenation, and possibly destruction, in NGC 7023 at radii 
< 30" (~ 0.058 pc) from HD 200775. The onset of PAH dehydrogenation in NGC 7023 
coincides with a drop in PAH intensity in every band and an ionization fraction of /+ ~ 
0.74. At projected distances > 30" the dehydrogenation ratio reaches a minimum before 
increasing again in the H2/PAH bright ribbon feature. The dehydrogenation signature in 
the molecular layer may indicate that H 2 is forming on the PAH surfaces, and that the 
process is removing a bound hydrogen atom from the PAH. We find that in all three PDRs, 
^12.7/^11.3 peaks where T gas ~ 600-700K, possibly indicating the peak efficiency temperature 
for H 2 formation on PAH grain surfaces. 

In three of the four PDRs studied, the PAH emission peaks a few arcseconds closer 
to the central star than the H 2 emission. The projected separation distance varied from ~ 
0.6" and 1.5" in NGC 2023S and NGC 7023 respectively, to ~ 10" in NGC 2023N. In IC 
63, no separation was observed between H 2 and PAH features, however the weak G /hh of 
IC 63 or the viewing geometry may make the separation unresolvable. In NGC 2023N we 
observed no significant separation between the PAHs and the ortho-H 2 , however the para-H 2 
dominates the total H 2 emission and column density in this PDR. This may also be the case 
in NGC 7023 where an ortho- and para-H 2 separation was observed of similar magnitude 
(~ 2") to the PAH-H 2 total separation. We lack the resolution to definitively determine if 
the ortho-H 2 is cospatial with the peak PAH emission in this PDR. No separation between 
ortho- and para-H 2 was detected in NGC 2023S, although there is a suggestion of a PAH-H 2 
separation of ~ 0.6", although it is not statistically significant. Higher spatial resolution 
measurements are required to confirm a separation. 
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6.2. H 2 Features 

Despite the strong PAH features in these PDRs, we were able to construct spectral line 
maps of multiple H 2 emission lines. We fit gas temperatures of the warm H 2 component 
which closely match global H 2 derived temperatures from previous studies, but with the 
added benefit of spatial resolution. Due to the strength of the PAH features, we were 
unable to accurately fit the H 2 S(4) and S(6) rotational lines in NGC 2023 and NGC 7023, 
however we estimate that the error introduced to the S(4) flux is negligible in IC 63. For 
the purposes of temperature and column density calculations we applied a correction to the 
S(4) flux to increase our spatial coverage area. 

In the bright emission ridges of NGC 7023 and NGC 2023N the ortho-rotational 
emission peaks before the para- rotational H 2 ; an effect most pronounced in NGC 2023N 
and not detected in IC63 and NGC 2023S. In NGC 2023N, where we had complete 
complementary IRS LL data, the peak of the ortho-to-para ratio calculated for the cool, 
low density H 2 , defined by the S(0) and S(l) emission lines, is cospatial with the peak of 
the ortho-to-para ratio for the warm, dense clumps; suggesting that H 2 self-shielding occurs 
at the same location regardless of the density or temperature phase of the H 2 gas. 

We presented evidence of the preferential self- shielding of ortho-H 2 in both NGC 7023 
and NGC 2023N. We also detected an implied trend in the ortho-to-para ratio in IC 63, with 
near LTE values at the leading edge of the PDR falling to ~ 1.5 at the edge of our coverage 
area. In NGC 7023 we measured Rop > 3 (Figure 0] (g)), however the measurements are 
not statistically significant and are consistent with the LTE value of Rop = 3. Ortho- and 
para-H 2 emission profiles in NGC 2023N and NGC 7023 (Figures [TT] and [TUT) show that in 
both PDRs the ortho-H 2 emission peaks starwards of the para-H 2 emission profile, which 
supports a preferential self- shielding scenario. In NGC 2023N, we observe no lasting effect 
of preferential self-shielding of ortho-H 2 on Rqp, as Rqp was nearly identical on either side 
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of the dissociation front. The ortho-to-para ratio of the cool H 2 component in NGC 2023N 
is consistent with an initial formation temperature of 50 - 75K. 

The authors would like to thank Dr. David Neufeld for helpful discussions and advice. 
We would also like to thank the anonymous referee for their helpful comments. This work is 
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Table 1. Literature Properties of Sample PDRs 



NGC 7023 NGC 2023N (pos.l) NGC 2023S (pos.2) IC 63 



Central Star 


HD 200775 


HD 37903 


HD 37903 


7-Cas 


Spectral Type 


B2.5 Ve a 


B1.5 V b 


B1.5 V b 


BO. 5 IVe c 


Te// (K) 


17000 a 


22000 b 


22000 b 


25000 c 


TVs (K) 


400-800 d 


550 


550 


620 c 


n H (cm -3 ) 


10 3 -10 6d 


10 4 -10 6c 


10 3 -10 6e 


10 4 -10 5c 


Go 


2400 d 


2200 c 


3-4000 c 


650 c 



Rogers et al 



(1995) 



Hde Boerl |l983j) 



J arisen et al 



I199J); 



Karr et al 



(2005) 



c Fuente et al 



(2000) 



Lemaire et al 



(1996); 



Chokshi et al 



(1988) 



Draine fc Bertoldi 



(1996); 



Steiman-Cameron et al 



(1997) 
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Table 2. Summary of Spitzer Observations 



Object 


Obs. 


Program 


Target 


Target 


Exposure 


Number of 


Background 


Background 




Date 


ID 


RA 


DEC 


Time (s) 


exposures 


Type 


ID 


NGC 7023 


Aug. 8, 2004 


28 


21:01:31 


+68:10:43 


14.68 


30 


Nearby Obs. 


22 


NGC 2023N 


Mar. 12, 2005 


3512 


05:41:35 


-02:12:41 


14.68 


21 


Nearby Obs. 


172 


NGC 2023S 


Oct. 8, 2007 


30295 


05:41:37 


-02:16:45 


29.36 


108 


SKY 




IC 63 


Jan. 7, 2005 


3512 


00:59:00 


+60:53:10 


14.68 


31 


Outrigger 
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Table 3. Average Properties for Regions with Calculated T, 



Object 


Distance 


T avg 


N(H 2 ) a „ 9 


Avg. 


Avg. 


Ortho-para 




from Star (") a 


(K) 


(10 20 cm- 2 ) 


I7.7/I11.3 


G /n H 


Ratio 


NGC 7023 


41.9 ± 11.7 


873 ± 273 


0.54 ± 0.53 


6.06 ± 1.74 


0.39 ± 0.38 


1.96 ± 0.67 


NGC 2023N 


174.8 ± 14.4 


522 ± 35 


0.94 ± 0.33 


5.53 ± 0.84 


0.32 ± 0.12 


1.64 ± 0.54 


NGC 2023S 


67.0 ± 24.0 


689 ± 149 


0.58 ± 0.40 


5.38 ± 1.00 


0.27 ± 0.10 


1.26 ± 0.57 


IC 63 


1192 ± 8.5 


637 ± 51 


0.26 ± 0.10 


3.23 ± 0.38 


0.10 ± 0.02 


1.99 ± 0.37 



Note. — ± values represent standard deviations of the values calculated, not the error 
a Projcctcd distance 
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Arcseconds 

Fig. 1— IRAC 3.6/mi, 8/im, and MIPS 24/im image of NGC7023 with IRS slit overlays. 
Overlapping slits are merged. North is oriented up in this image. 
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NGC2023 
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Fig. 2.— IRAC 3.6/mi, 8/im, and MIPS 24/zm image of NGC2023 with IRS slit overlays. 
Overlapping slits are merged. North is oriented up in this image. 
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Fig. 3. — IRAC 3.6/im, 8/im, and MIPS 24/im image of IC63 with IRS slit overlays. The 
star 7 Cas is approximately 1.3 pc SW of this emission nebula and not in this image. North 
is oriented up. 
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Fig. 4. — PAH and H2 emission and diagnostics for NGC 7023. (a) Total detected PAH 
emission, (b) Total detected H 2 emission, (c) /7.7//11.3. (d) 7 An. 3- ( e ) H2 column 
denisty. (f) H 2 derived gas temperature, (g) H 2 Ortho-to-para ratio, (h) G /uh- 



NGC 20$3 North 



PAH Emission 




H 2 Emission 




-20-10 10 20 



-20-10 10 20 



( c ) - !7.7/l 





0.60 , 



-20-10 10 20 



-20-10 10 20 



(e) " N(H 2 ) 




(f) " T 




-20-10 10 20 



-20-10 10 20 



(g) - Rq 




(h), -G /n H 




-20-10 10 20 



-20-10 10 20 



Fig. 5. — PAH and H2 emission and diagnostics for NGC 2023 North, (a) Total detected 
PAH emission, (b) Total detected H 2 emission, (c) /7.7//11.3. (d) h2.7A11.3- ( e ) H 2 column 
denisty. (f) H 2 derived gas temperature, (g) H 2 Ortho-to-para ratio, (h) G /uh- 
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Fig. 6. — PAH and H2 emission and diagnostics for NGC 2023 South, (a) Total detected 
PAH emission, (b) Total detected H 2 emission, (c) /7.7//11.3. (d) Z12.7A11.3- (e) H 2 column 
denisty. (f) H 2 derived gas temperature, (g) H 2 Ortho-to-para ratio, (h) G /uh- 
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Fig. 7. — PAH and H2 emission and diagnostics for IC 63. (a) Total detected PAH emission, 
(b) Total detected H 2 emission, (c) /7.7//11.3. (d) /12. 7/^11.3- (e) H 2 column denisty. (f) H 2 
derived gas temperature, (g) H 2 Ortho-to-para ratio, (h) G /uh- 
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Fig. 8. — Plot of the H 2 S(0) - S(5) rotational lines as a function of transition energy 
for a single data cube element near the bright ribbon of NGC 7023. Overplotted are two 
temperature fits, one for a higher column density gas (3.5 x 10 22 cm~ 2 ) at 125K (dot-dash) 
and another for a lower column density gas (9 x 10 19 cm" 2 ) at 650K (dash). 
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Fig. 9. — Average normalized emission profile for PAHs (solid), H2 (dash) and (I7.7/In.3)/12 
(dot) for NGC 7023. Multiple cross-sections through the PDR originating from the direction 
of HD 200775 are averaged to produce this profile. 
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Arcseconds 

Fig. 10. — Plot of the ortho-H 2 S(3) (dash) and para-H 2 S(2) (solid) emission in NGC 7023 as 
a function of projected distance from the central star. Also plotted is the average normalized 
temperature profile (dash-dot) and the ortho-para ratio (dot). 



-48- 




Fig. 11. — Plot of the total ortho-H 2 (dash) and para-H 2 (solid) emission in NGC 2023N as 
a function of projected distance from the lower right corner of the maps in Figure |5j Also 
plotted is the average normalized temperature profile (dash-dot) and the ortho-para ratio 
(dot). 
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Fig. 12. — Average normalized emission profile for PAHs (solid), H2 (dash) and (I7.7/In.3)/12 
(dot) for NGC 2023S. Multiple cross-sections through the PDR originating from the direction 
of HD 37903 are averaged to produce this profile. 
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Fig. 13. — Average normalized emission profile for PAHs (solid), H2 (dash) and (I7.7/In.3)/12 
(dot) for IC 63. Multiple cross-sections through the PDR originating from the direction of 
7-Cas are averaged to produce this profile. 
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Fig. 14. — Plot of the average profiles of I12.7/A1.3 (triangle) and N(H 2 ) (asterix) versus 
distance to HD 200775 in NGC 7023. 



